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c/3 ■ Some relationships, fundamental to the resolution of interface wave 



problems, are presented. These equations allow for the derivation of 
explicit secular equations for problems involving waves localized near 
the plane boundary of anisotropic elastic half-spaces, such as Rayleigh, 

'^ . Sholte, or Stoneley waves. They are obtained rapidly, without recourse 

to the Stroh formalism. As an application, the problems of Stoneley 
O ' wave propagation and of interface stability for misaligned predeformed 

incompressible half-spaces are treated. The upper and lower half- 
spaces are made of the same material, subject to the same prestress, 
*zl I and are rigidly bonded along a common principal plane. The principal 

Q>^ ' axes in this plane do not however coincide, and the wave propagation 

is studied in the direction of the bisectrix of the angle between a 

^~ I principal axis of the upper half-space and a principal axis of the lower 

^ ■ half-space. 
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1 Introduction 

We all know from experience or from intuition that when we press together 
two large, well-polished, glass plates, they will stick together extremely well; 
in fact we will have a tough job trying to separate them again because in 
effect, the two glass panes have become one. A somewhat similar process 
of "gluing without glue" is used in the microelectronics and optoelectronics 
industries to bond together semiconductor wafers (Gosele and Tong, 1998). 
When brought into contact, mirror-polished, flat, clean wafers made of al- 
most any material are attracted via Van der Walls forces and adhere in a 
rigid and permanent way. This method of direct bonding allows for new and 
promising designs for insulators, sensors, actuators, nonlinear optics, light- 
emitting diodes, etc. Solid polymers can also be brought into permanent and 
rigid contact to manufacture polymer composites. The main traditional tech- 
nologies of polymer joining are: mechanical fastening (bolts, rivets, fit joints) 
and: adhesive bonding. Another technology, "fusion bonding" , presents great 
advantages over the previous ones such as, avoidance of high stress concen- 
trations, reduced surface treatment, less inhomogeneities at the interfaces, 
etc. A recent book by Ageorges and Ye (2002) presents a comprehensive 
description of fusion bonding, defined as "the joining of two polymer parts 
by the fusion and consolidation of their interface" ; in particular four classes 
of fusion bonding are listed: "bulk heating (co-consolidation, hot-melt adhe- 
sives, dual- resin bonding), frictional heating (spin welding, vibration welding, 
ultrasonic welding), electromagnetic heating (induction welding, microwave 
heating, dielectric heating, resistance welding), and two-stage techniques (hot 
plate welding, hot gas welding, radiant welding)." 

For semiconductor wafer bonding, the most common and most economi- 
cal combination is the silicon/silicon wafer bonding. Mozhaev et al. (1998) 
considered the theoretical implications of misorientation when two identi- 
cal silicon wafers are rigidly bonded. Specifically they considered, within the 
framework of anisotropic linear elasticity, the propagation of interface (Stone- 
ley) waves along the bisectrix of a twist angle of misorientation . This paper 
deals with the propagation of waves of similar nature but within the frame- 
work of small motions superimposed on large static deformations (stress- 
induced anisotropy) in nonlinear elasticity, which is required to describe the 
possible large elastic deformations of polymers (elastomers). Specifically, two 
half-spaces made of the same hyperelastic incompressible material are main- 
tained in the same static state of pure homogeneous deformation and are then 
rigidly bonded along a common principal plane of deformation, but in such a 
way that the two principal axes defining this plane for one half-space do not 
coincide with the two counterpart principal axes for the other half-space. 



The study of a superimposed infinitesimal interface wave propagating 
along the bisectrix of the misalignement angle provides insights into the 
possible ultrasonic non-destructive evaluation of the bond and of the an- 
gle of twist, and into the influence of the pre-strain on the interface stability. 
Indeed in general, an increasingly tensile load applied on the semi-infinite 
bodies leads to faster speeds for an interfacial wave whereas a compressive 
load slows the wave down, until eventually the buckling/bifurcation crite- 
rion is met at the critical load, where the speed is zero. Such investigations 
of interface waves and interfacial stability for two bonded hyperelastic pre- 
strained half-spaces are quite rare in the finite elasticity literature. They 
were initiated by Biot (1963), followed by a handful of papers: Chadwick & 
Jarvis (1979a, 6), Dunwoody & Villaggio (1988), Dowaikh & Ogden (1991), 
Chadwick (1995), Destrade (to appear, a). Note however that in these ar- 
ticles, the two semi-infinite deformed bodies are always aligned so that all 
three principal axes of pre-deformation for the upper half-space coincide with 
those for the lower half-space. Moreover, the interface wave propagates along 
one of the common principal axes, except for Chadwick & Jarvis (1979a, 6), 
who consider non-principal directions of propagation but for a specific simple 
form of strain energy function (compressible neo-Hookean materials). 

Some papers have been devoted to the study of non-principal interface 
waves (e.g. Flavin, 1963; Willson, 1973, 1974; Chadwick & Jarvis, 1979a, b, c; 
Connor & Ogden, 1995; Rogerson & Sandiford, 1999; etc.) but explicit sec- 
ular equations for tri-axially pre-strained materials have been found only in 
the case of neo-Hookean (compressible or incompressible) materials. This 
strain energy density is exceptional with respect to the propagation of inho- 
mogeneous plane waves because for any direction of propagation in a princi- 
pal plane (with associated orthogonal attenuation), the in-plane strain com- 
ponents always decouple from the anti-plane strain components. Section 
2 presents the incremental equations of motion for non-principal interface 
waves in a tri-axially pre-strained material with a generic (non neo-Hookean) 
strain energy function, formulated as a first-order differential system for the 
six-component displacement-traction vector. It also presents the "effective 
boundary conditions" , consequences of the continuity of the mechanical dis- 
placements and tractions across the interface. Namely, these conditions are 
that at the boundary, either one displacement and two tractions are zero, or 
two displacements and one traction are zero. Mozhaev et al. (1998) noted 
this important result for bonded silicon/silicon wafers in linear anisotropic 
crystallography. It is further proved here that once the displacement-traction 
vector is normalized with respect to one of its non- zero components, then 
the normalized components are either real or pure imaginary (see Appendix 
for details). This Author recently obtained, in a quick and simple manner. 
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Figure 1: Misoriented deformed cuboids at the interface 



some equations which are fundamental to the resolution of general interface 
boundary value problems (Destrade, 2003; Destrade, to appear, b) includ- 
ing Rayleigh (solid/vacuum interface), Scholte (solid/fluid interface), and 
Stoneley (solid/solid interface) waves; they are presented in Subsection 3.1. 
The next Subsection brings together these fundamental equations and the 
effective boundary conditions for the title problem, with an explicit form of 
the secular equation that is, a polynomial of which the interface wave speed 
is a root. Finally Section 4 shows how numerical results can be obtained 
from the analysis, with the case of deformed semi-infinite bodies made of 
Mooney-Rivlin material, a model often used to describe the behaviour of 
incompressible rubber in large deformations. 



2 Basic equations 

2.1 Finite static pre- deformation 

Consider an infinite body made of two pre-strained liyperelastic semi-infinite 
bodies, rigidly bonded along a plane interface. Take the origin O of a rectan- 
gular Cartesian coordinate system (OXYZ) to lie in the boundary, so that 
the plane of separation between the half-spaces is F = 0. The upper half- 
space y > is made of an incompressible isotropic hyperelastic body, with 
mass density p and strain energy function W, which has been subjected to a 
finite static pure homogeneous deformation with principal stretches Ai, A2, A3 
(Ai 7^ A2 7^ A3 7^ Ai and A1A2A3 = 1) say, along the principal axes OX1X2X3 
say, such that X2 is aligned with Y but Xi (X3) makes an angle 6 with X 
(Z). The lower half-space F < is made of the same body, which has been 
subjected to the same pre-deformation that is, a pure homogeneous defor- 
mation with principal stretches Ai, A2, A3 along the principal axes OX1X2X3 
(say) where X2 is aligned with Y, but such that now Xi (X3) makes an angle 
—6 with X {Z). Figure 1 summarizes this set-up with the representation of 
two parallelepipeds, one above Y = 0, one below, which were unit cubes 
before the static deformation took place. 

2.2 Incremental equations of motion and effective bound- 
ary conditions 

Now a small-amplitude inhomogeneous plane wave is superimposed upon 
the primary large static deformation. The wave propagates with speed v and 
wave number k in the X direction and vanishes away from the interface Y = 
0. In other words, the corresponding incremental mechanical displacement 
u is of the form, 

u(x, y, z, t) = U(%)e^'=(^-^*), U(±oo) = 0, (1) 

where [Oxyz) = (OXYZ) is the rectangular Cartesian coordinate system 
associated with the motion. 

Similarly, the incremental nominal tractions 5^2 = 1? 2, 3) acting on the 
planes y = const, are of the form, 

s,2{x, y, z, t) = ikt,{ky)e"'^^-'"\ t,(±oo) = 0. (2) 

The following quantities allow for a compact form of the incremental 



equations of motion (Rogerson and Sandiford, 1999), 

2A, := \lW,, - 2\\W,, + \]Wj, + 2{\W, - \W,)\\/{\^, - A^^) = 2/3,,, 

(3) 



(where Wi := dW/d\i) and also, 

C0 := cos^, S0 

r]0 := 2cg(/3i2 + 721) + Se73i, t^e 
l^e ■= cl-fn + 2sj{l323 + 723), K,e 



sin6, 

£9(712 - 721) + Se(732 - 723), 

C0Sg{f3i3- P12- f323- 121-123)- (4) 



Then a careful reading of the literature on incremental motions in in- 
compressible materials (e.g. Ogden, 1984; Chadwick et ai, 1985; Chadwick, 
1995; Rogerson & Sandiford, 1999; Destrade & Scott, to appear) and a fair 
amount of algebra reveal that the incremental equations of motion (in the 
upper half-space) can be cast in the form. 



^' = iN|, where ^{ky) = [U{ky),t{ky)f 



(5) 



the prime denotes differentiation with respect to ky, and N is the following 
matrix, 



N 







-C0 




-C0 








-S0 


i 

721 











X 




-V9 




-S0 








K0 








-C0 


1 

723 









X 


-P0 







-C0 





-S0 




K0 







X 


-^e 





-S0 






X ■-- 



pv"^. 



(6) 



The block structure of this matrix is reminescent of that for the "fundamental 
elasticity matrix" N of linear anistropic elasticity (Ingebrigsten and Tonning, 
1969), although its components are different. In the lower half-space the 
equations of motion are the same, allowing for the change in sign of S0 and 

K0. 

At the interface, the displacement-traction vector is continuous. 



|(0+)=^(0-)=:|(0). 



(7) 



Lengthy, but straightforward to establish, expressions for these quantities 
lead to the effective boundary conditions at the interface. In effect, and as 
displayed by Mozhaev et al (1998) in the linear anisotropic elasticity case. 
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either one displacement component and two tractions are zero at the inter- 
face, or vice-versa. More specifically, and these details were not noted by 
Mozhaev et al (1998), ^(0) must be of one of the two following forms. Either 

|(0) = [/2(0)[0,l,m2,/3i,0,0]T, (8) 

where 0^2 and /3i are real, or 

^(0) = t/i(0)[l,0,0,0,«i,^/32]^, (9) 

where ai and (32 are real. The Appendix presents the derivation of these 
expressions. Following Mozhaev et al (1998), the interface acoustic wave 
satisfying ^ (resp. ^) is called lAWl (resp. IAW2). 

3 Explicit resolution 

3.1 Fundamental equations for interface boundary value 
problems 

Here a quick derivation is made of some equations (Destrade, 2003) which 
are fundamental to the resolution of interface boundary problems involving 
waves or static deformations which are localized near, or equivalently vanish 
away from, the plane interface of a semi-infinite body. The equations are 
valid for unconstrained or constrained, pre-deformed nonlinearly elastic or 
anisotropic linearly elastic, materials. They do not rely on the Stroh (1958) 
formalism. All that is required for their derivation are the few manipulations, 
available in the literature, leading to the equations of motion and boundary 
conditions in the form ([5]), ^2, ([2]) 2, that is 

^' = zNt |(oo) = 0. (10) 

In general here, N is a square matrix of even dimensions, 2p x 2p say. 
Now a simple induction process (Currie, 1979; Ting, 2003) shows that N", 
where n is any positive or negative integer, is of the following form (see (jS]) 
for n = 1), 



^n 



K(°) Ni(^)T 



with K(°)=K(-)T, Na^") = Na^")^. (11) 



Consequently, pre-multiplication of N" by I, defined as 



I 



Ip 
Ip 



where Ip is the p x p identity matrix, (12) 



leads to a symmetric matrix, 



IN^ 






(IN 



n\T 



(13) 



Now take the scalar product of flTU]) i by IN*^^ and add the complex con- 
jugate quantity to obtain ^'-IN^^ + ^-IN"^ = 0. Direct integration between 
zero and infinity yields, using (fTOjl o. the fundamental equations sought. 



^(0)-iN"^(0) = 0. 



(14) 



The Cayley-Hamilton theorem states that N satisfies its own characteristic 
polynomial of degree 2p. Hence there are only 2p — 1 powers of N which 
are linearly independent and so, Eqs. (IT^ generate at most 2p — 1 linearly 
independent equations. For instance in our case N is a 6 x 6 matrix, leading 
to only 5 linearly independent equations. These are nevertheless enough to 
deduce the secular equation explicitly, as is now seen. 



3.2 Secular equation 

The fundamental equations flT4|) now prove useful in the resolution of the 
paper's problem. The structure of IN" depends on the parity of n. For 
instance. 
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for n 



-2,2, 



(15) 



where NJ'j := (N")^^; thus for example, Nj^ := {1^^)42 = iV'4fciVfc2- The forms 
([8]) and ([9]) for ^(0) lead to two trivial identities when flTH) are written at 
n = —2, 2. On the other hand. 
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for n = —1, 1, 3, 



(16) 



and so, the fundamental equations (IT^ . written for lAWl that is, ^(0) given 
by dH]), and for n = —1, 1, 3, yield the following system. 
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(17) 



When the fundamental equations fll4p are written for IAW2 that is, ^(0) 
given by (jH]), and for n = —1,1,3, they yield the system. 
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These equations are of the type Mi^bk = ca where here M is the 3x3 
matrix on the left hand-side of ([IIj) for lAWl and of ^) for 1AW2, a is 
the vector on the right hand-side, and b = [2/3i,/3^, af]'^, [2a;i, af,/3|]'^, for 
lAWl and for IAW2, respectively. By Cramer's rule, the unique solution is 
hk = Afc/A, where A = det M, and A^ is the determinant of the matrix 
obtained from M by replacing the k-th vector column with a. However 
two components of b are related one to another. Specifically, b1 = 4^2, 
which is the explicit secular equation for interface waves propagating along 
the bissectrix of the mis orientation angle for two identical, rigidly bonded, 
pre- deformed, hyperelastic half-spaces, 

A2-4AA2 = 0. (19) 



4 Numerical results for rubberlike materials 

The classical Mooney-Rivlin model has been used extensively to model the 
behaviour of incompressible rubberlike materials in large deformations. Its 
strain energy function W is given by 

2W = C{Xl + A^ + A^ - 3) + D{XlXl + XJXl + A^A? - 3), (20) 

where C and D are constant material parameters. In that case, the 7jj and 
f3ij defined in (^ reduce to {i ^ j ^ k ^ i), 

7., = A,2(C + DXl), 2A, = {X} + X]){C + DXl) = 7., + 7,.. (21) 

Now the peculiarities of interface acoustic waves in deformed hyperelastic 
materials are highlighted. In particular, and in contrast with the correspond- 
ing situation in linear anisotropic elasticity, it is seen that there exist certain 
ranges of misorientation and certain ranges of pre-stretch ratios for the exis- 
tence of an lAW. Also, interface instability may arise at certain compressive 
critical ratios. 
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4.1 Mooney-Rivlin material in tri-axial strain 

First, in order to make the connection with results obtained by Rogerson & 
Sandiford (1999) about non-principal surface waves, the material constants 
C, D, and the stretch ratios A^ are fixed at the following values, 

C = 2.0, D = 0.8, \l = 3.695, Xj = 0.7, A^ = 0.387. (22) 

The secular equation ( 1T9|) is a polynomial of degree 10 in X = pv'^ for 
lAWl, of degree 6 for IAW2. Out of the 16 possible roots, only one leads to an 
interface wave of the Stoneley type satisfying the known conditions that its 
speed is bounded above by the speed of the slowest homogeneous bulk wave 
and below by the speed of the Rayleigh surface wave associated with either 
half-space (Barnett et ai, 1985). Moreover this wave, of the lAWl type, 
satisfies these requirements only within a limited range of misorientation 
angle, approximatively 16.7° < 6 < 73.6°. This situation is in sharp contrast 
with the case of silicon/silicon wafers in linear anisotropic elasticity (Mozhaev 
et ai, 1998) where the lAWl was found to exist for all 6. 

Figure 2 depicts the variations of the relevant root to the secular equation, 
scaled as vX = a//w^, with 9 (thick curve). The speeds of two homogeneous 
shear waves (roots of det N = 0) are represented by the two thin curves 
above, and the speed of the Rayleigh wave, by the thin curve below. 

4.2 Mooney-Rivlin material in compressive plane strain 

Next, the half-spaces are assumed to have been pre-deformed in such a way 
that they were not allowed to expand in the Z-direction (A3 = 1). For this 
example, the different parameters take the following numerical and algebraic 
values, 

^ = 2.0, D = 0.8, = 30°, Ai = A, A2 = A"\ A3 = 1. (23) 

Figure 3 shows the variations of the speeds as functions of A in compres- 
sion (A < 1). The two upper dashed curves represent the speeds of homoge- 
neous shear waves; the intermediate thick curve, the speed of the Stoneley 
wave (type: lAWl); the bottom thin curve, the speed of the Rayleigh wave 
associated with either deformed half-space. Again, the situation is differ- 
ent from that encountered in linear anisotropic elasticity with silicon/silicon 
wafers (Mozhaev et ai, 1998). In particular, the Stoneley wave exists only 
for stretch ratios greater than 0.321 and the Rayleigh wave for stretch ratios 
greater than 0.451. In between these two critical stretches, there is a range 
were the Stoneley wave exists but not the Rayleigh wave. At the critical 
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Figure 2: Mooney-Rivlin bimaterial subject to the triaxial pre-stretch A^ = 
3.695, A2 = 0.7, A| = 0.387: plot of the Stoneley wave speed with 6 (thick 
curve), bonded above by the shear wave speeds (crossing at ^ = 45°) and 
below by the Rayleigh wave speed. 
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Figure 3: Deformed Mooney-Rivlin bimaterial in compressive plane strain 
at 6 = 30°: variations with A of the shear (dashed curves), Stoneley (thick 
curve), and Rayleigh (lower curve) wave speeds. 

stretch, instability might occur; hence is presented an example where two 
bonded deformed half-spaces exhibit interfacial stability at some compres- 
sive stretch ratios for which the separated half-spaces are unstable (at least 
in the linearized theory). 
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Appendix: effective boundary conditions 

Here two main results about the boundary conditions of the paper's problem 
are established. The first result is that at the interface, either one displace- 
ment component and two tractions are zero, or vice-versa. Mozhaev et al. 
(1998) found this result for the corresponding problem in anisotropic linear 
elasticity; because of space limitations they only sketched the proof. The 
second result is that once the displacement-traction vector at the interface is 
normalized with respect to one of the three non-zero components, then the 
two others are either purely real or purely imaginary quantities. This Ap- 
pendix is a generalization of results obtained by Ting (to appear) for surface 
(Rayleigh) waves and by Destrade (2003) for 2-partial interface (Stoneley) 
waves. 
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Consider, for the upper half-space, a combination of inhomogeneous plane 
waves of the form, 

where ^* is a constant vector and '^{pi) > to ensure decay of the wave with 
distance from the interface. Substitution of this form of solution into the 
equations of motion ([5]) shows that the pi are roots of the bicubic: det (N — 
pi) = 0, or (Rogerson and Sandiford, 1999), 

721723/- [(721+723)X-Cl]/+(X2-C2X + C3)p' + (X-C4)(X-C5) =0, (2) 

with 

Cl = (721713 + 2/3i2723)Ce + (723731 + 2/32372l)Se, 
C2 = (723 + 7l3 + 2/3i2)Ce + (721 + 731 + 2/323) -5^, 
C3 = (712723 + 2/3i27l3)c0 + (721732 + 2/32373i)s0 

+ [712721 + 713731 + 723732 - (/3l3 " /3l2 " /323)^ + '^I^12l^2i]clsl, 
Ca = 7l2Ce + 732S^, 

C5 = 713C0 + 2(3i3clsg + 73i4- (3) 

Let pi,P2,P3 be the three roots with positive imaginary part. Then the 
wave combination is 

£ = ry^il^iK^+Piy--"*:) _|_ ry^t2^ik{x+p2y-vt) _|_ ^ tS^ik{x+p3y-vt) / ^ g\ /^\ 

where the 7^ are constant numbers. 

In the lower half-space, the equations of motion are of the same form as 
in the upper half-space with 6 changed to —6. Hence only sq and k change 
signs in the components of N whilst the coefficients of the associated charac- 
teristic polynomial remain unchanged. Consequently, a suitable combination 
of inhomogeneous plane waves is here 

I = ^^|ie^'=(^-Pi2'-''*) + ^2|2e*'=(^-P2J^-^*) + ^3|3gifc(x-p3s/-i't)^ (y < 0), (5) 

where the 7^ are constant numbers and the ^* constant vectors. 
Now two possibilities arise for the roots Pi,P2,P3 of ([2]). Either 

(I) Pi = iqi {qi > 0), or (II) pi = -p^, ps = iq^ {q^ > 0). (6) 

In Case (I) it is a straightforward matter to show, using the adjoint of 
(N — pi), that the ^* (upper half-space) can be written in the form, 

C = [ai, ihi, Ci, idi, Ci, ifif, (7) 
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say, where ai,bi,Ci,di,ei, fi {i = 1,2,3) are real numbers. Furthermore, the 
vectors ^* (lower half-space, where 6 is changed to —6) are then in the form, 

C = [-fli, ibi, Ci, idi, -Ci, -ifi]^. (8) 

At the interface y = the displacement-traction vector is continuous, $,{x,0, t) 
^(a;,0,t), or 

71^' + 72^^ + 73^^ = 7i|' + l2i'' + 73!^ 



These six equations are recast as 

A(7 + 7) = 0, B(7-7) 
where 



0, 



(9) 
(10) 



A 



ai 


02 


03" 


ei 


62 
^/2 


63 



B 



ihi 


ib2 


ib'i 


C\ 


C2 


C3 


id\ 


id2 


id^_ 



1 



"71" 




71" 


72 


, 7 = 


72 


73. 




73. 



11) 



For non-trivial solutions to exist, either (a) det A = or (b) det B = 0. In 
case (a), Eq. (ITOJI q leads to 7 = 7. Then Eq. (ITOl) i reads 



A7 = 0. 



(12) 



Owing to the form flTT|) i of A, this condition is satisfied when the components 
of 7 are all real. In conclusion, ^(0) now reads 



m = 71 



"0" 




"0" 




"0" 


ihi 




ih2 




ib3 


%d\ 


+ 72 


C2 
id2 


+ 73 


C3 

ids 















_ 0_ 




_0_ 




_ 0_ 



7i real. 



(13) 



that is, ^(0) is of the form ([8]). A similar procedure shows that in case (b), 
^(0) is of the form Q. 

In Case (II) the ^* (upper half-space) can be written in the form. 
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ai 




en 




as 


bi 




-fei 




ibs 


Cl 

d. 


, e = 


Cl 


, e = 


C3 
ids 


ei 




el 
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/i. 




-/i. 




}h 



(14) 



16 



say, where ai,bi,ci,di,ei, /i are complex and 03, 63, C3, cis, 63, /3 are real. The 
vectors ^* (lower half-space) are then in the form, 



i' 



bi 

Cl 

di 
-ei 
-/i 



— «! 

— (ii 
. /i. 



-as 

ibs 

C3 

2^3 

-63 



(15) 



The continuity of the displacement-traction vector at y 
written in the form fITT]) where now 



can again be 
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ai 
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as 


ei 


el 
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/i 


-/i 


^/s 
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bi 
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Cl 


cl 


cs 


rfi 


-rfl 


4. 



(16) 



Again case (a) det A = or case (b) det B = arise. In case (a), 7 = 7 
and so A7 = 0. Owing to the form (I16j) i of A, this condition is satisfied 
when 7 is of the form 7 = [71,7157s]"'") where 71 is complex and 73 real. In 
conclusion, ^(0) now reads 



m = 71 



"0" 




" " 




" 0" 


bi 




-bi 




ib^ 


Cl 

rfi 


+ 71 


-rfi 


+ 7s 


Cs 

id'i 

























_0_ 



7s real. 



(17) 



which means that, once normalized, ^(0) is of the form 
dure shows that in case (b), ^(0) is of the form ([9]). 



). A similar proce- 



17 



